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ABSTRACT 

We test the assumption of strict hydrostatic equilibrium in galaxy cluster MS2137. 3-2353 (MS 2137) 
using the latest CHANDRA X-ray observations and results from a combined strong and weak lensing 
analysis based on optical observations. We deproject the two-dimensional X-ray surface brightness 
and mass surface density maps assuming spherical and spheroidal dark matter distributions. We 
find a significant, 40%-50%, contribution from non-thermal pressure in the cor e assuming a spherica l 
model. This non-thermal pressure support is similar to what was found by iMolnar et all ([20101 ) 
using a sample of massive relaxed clusters drawn from high resolution cosmological simulations. We 
have studied hydrostatic equilibrium in MS 2137 under the assumption of elliptical cluster geometry 
adopting prolate models for the dark matter density distribution with different axis ratios. Our results 
suggest that the main effect of cllipticity (compared to spherical models) is to decrease the non-thermal 
pressure support required for equilibrium at all radii without changing the distribution qualitatively. 
We find that a prolate model with an axis ratio of 1.25 (axis in the line of sight over perpendicular to 
it) provides a physically acceptable model implying that MS 2137 is close to hydrostatic equilibrium 
at about 0.04-0.15 K V i r and have an about 25% contribution from non-thermal pressure at the center. 
Our results provide further evidence that there is a significant contribution from non-thermal pressure 
in the core region of even relaxed clusters, i.e., the assumption of hydrostatic equilibrium is not valid 
in this region, independently of the assumed shape of the cluster. 



Subject headings: galaxies: clusters: individual (MS2137.3-2353)-gravitational lensing: 
gravitational lensing: weak -X-rays: galaxies: clusters 
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1. INTRODUCTION 

Galaxy clusters are the largest gravitationally bound 
and virialized systems in our universe. They are im- 
portant cosmological tools to probe the curvature, mass 
and energy density of the universe. One of the most 
important parameters of galaxy clusters for cosmology is 
the total mass. Cosmological simulations make specific 
predictions for the total mass distribution of clusters 
and the distribution of mass in individual clusters. The 
total mass of clusters is a particularly important scaling 
parameter when using galaxy clusters to constrain 
cosmological parameters. Unfortunately, the mass of 
clusters can not be observed directly, it has to be derived 
from observations with some additional assumptions. 
The most widely used (and independent) methods to 
determine cluster masses are based on gravitational 
lensing and X-ray observations. The method based on 
gravitational lensing is independent of dynamics because 
the masses is derived from the tangential shear of the 
background galaxies distorted by the gravitational po- 
tential of the cluster. On the other hand, cluster masses 
estimated from X-ray observations is based on dynamics. 
The X-ray method assumes that the gas is relaxed and 
in strict hydrostatic equilibrium (i.e. the cluster is 
supported only by the thermal pressure). However, 
there is a long-standing problem that these two meth- 
ods so metimes result in sign i ficant l y different cluster 
mass es dWu et all 119981: 1 Ailed 119981: faicks fc Canizaresl 
20001 iHoekstral 120071: iPiffaretti k Valdarninil 12005 
Peng et all l2009t [R iemer-S0rensen et al.1 



Meneghetti et al.1 120101: iZhang et all 120101: lOkabe et al 



megl 

20101; IMorandi k Limousin et all 120111: iMorandi et al 
2011al) . 

The most popular method to derive the total 
mass of relaxed clusters from X-ray observations as- 
sumes that they are in strict hydrostatic equilib- 
rium and have spherical symmetry. Unfortunately, 
there are several physical processes which would break 
the strict hydrostatic equilibrium in clusters, for ex- 
ample: turbulent gas motion, cosmic rays, heat- 
ing, cooling, and mag netic fields (lEnsslin et all 11997 : 
Norman fe Brvadll999t iRasia et al.N2004t iSarazid 2004 : 
Rasia et al.1 120061: iNagai. Vikhlinin k Kravtsovl 12007: 



Pfrommer et alll2007l 120081: lPfrommerll2008t iGuo k Oh 
20081: iSkillman et all [200g iFang et all 120091 : iLau et al 



2011 



200S ; Ivazza et al. 120091: iLaganaet al ■ M2010t IWdaTnhn l 



iParrish et aljl2012f ) 

Numerical simulations suggest that turbulent mo- 
tions of the intracluster gas (ICG) caused by merg- 
ers and shocks provide a significant non-thermal pres- 
sure support which varies as a function of radius 
(iNorman fe Brvad 119991 iDolag et all 120051: IVazza et all 
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Lau et all 120091 : IVazza et ; al.l 120091 : IMolnar et al Il2010^ 

Based on smooth particle hydrodynamics (SPH) simu- 
lations, IVazza et all (|2006l) showed that the turbulence 
energy of gas particles scales approximately with ther- 
mal energy of clusters. In the outskirts of clusters, the 
ICG is also supported by the bulk g as motion and tur - 
bulcnce created by accretion shocks ([Burns et al.ll2010| ). 
The scaling relations between the thermal and turbulence 
energy in the different rad ii among the differe nt redshifts 
have been investigated by IVazza et all ([20111 ) . Their re- 
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suits suggest that the non-thermal energy increases with 
radius and it may be m ore than 20% o f the thermal en- 
ergy near the R V i r . iParrish et al.l (|2012f ) conducted 
magnctohydrodynamics simulations to study the turbu- 
lences deduced from the magnetothermal instability in 
the outskirts of the clusters, and suggested that non- 
thermal pressure can provide 5 - 30% support against 
gravity beyond i?soo- lLau et all (|2009f ) studied the effect 
of gas motion in the mass estimation of simulated clus- 
ters, and concluded that the non-thermal pressure due 
to residual gas motion at large radii, if not taken into 
account, would result an underestimate of the cluster 
mass of about 15% confirming previous results. Based 
on an analysis of massive cluster s draw n from cosmo- 
logical simulations, iMolnar et al. I (|2010l ) suggested that 
the non-thermal pressure can provide a significant, about 
30%, support in the central region, a minimum support 
at about 0.1-0.2 K V i r , and incr e ases w ith radius to 
about 35% at R mr . iFang et al.l (|2009f ) suggest that 
the pressure from rotation and streaming motions are 
significantly comparable to the pressure raised from ran- 
dom turbulent pressure out to certain radius range. The 
hydrostatic equilibrium assumption could result an un- 
derestimate of the mass of clusters due to t he neglected 
contr ibution from turbulent bulk motion (|Rasia et al.l 
l2006h . The mass estimation of the cluster could be 
fairly reconstructed if the gas motion pressure is taken 
into account in the hydrostatic e quilibrium equation 
(|Rasia et al.l 120041 : lLau et al.ll2009| ). In addition, inap- 
propriate assumption of the gas and the dark matter halo 
distribut ions could also bias the mass estimation of the 
clust e rs dGavazz l27MICor¥ss fe King|l2008HOguri et al 



20101; IMorandi fe Limousin et all 120111: IMorandi et al' 



2?niallbl: ISereno fc Umetsul 120 lit ISereno et all I20ia 
Sereno fc Zitrinll2012l ). " 



On the other hand, gravitational lensing of background 
galaxies is an unique, direct probe of the distribution of 
matter in galaxy clusters, determining the mass profiles 
from two-dimensional (2D) projected physical quantities 
such as reduced shear for example. However, it is ex- 
tremely hard to fully reconstruct the three-dimensional 
(3D) mass distribution even with triaxial models fitting 
to the lensing data, because the lensing technique is sen- 
sitive to all the projec ted mass in the line of sight (LoS) 
(|Corless fc Kind |2008|) . Methods based on lensing as- 
suming spherical symmetry would possibly overestimate 
(underestimate) the cluster masses if the mass distribu- 
tions are prolate (oblate) with the non-de generated prin - 
ciple axis aligned in LoS (see for example lGavazzill2005D . 
Moreover, in order to derive accurate mass profiles at all 
radii, it is necessary to combine strong and weak lensing 
observations which are sensitive to small and large radii. 

Molnar et al. (2010, hereafter M10) studied Abell 1689 
using the latest CHANDRA data and three massive clus- 
ters of galaxies drawn from cosmological adaptive mesh 
refinement (AMR) simulations. Assuming spherical ge- 
ometry, they found a significant, about 45% and 35%, 
contribution from non-thermal pressure in the core re- 
gions of Abell 1689 and simulated clusters respectively. 
IMorandi fc Limousin et all (|2011| ) fitted triaxial dark 
matter mass models with their major axis aligned with 
the LoS to the full lensing and X-ray data simultaneously, 
and derived the dark matter axis ratio is 2.02 ± 0.01 (1.24 
± 0.13) along the LoS (in the plane of sky) for Abell 1689. 



But, unlike M10, IMorandi fc Limousin et al.l (|2011ft as- 
sumed that the non-thermal pressure contribution is con- 
stant throughout the all radii. They found that the non- 
thermal pressure contributes about 20% to the total pres- 
sure needed if th e hydrostatic equilibrium h o lds. U sing 
similar methods. IMorandi fc Limousin et"al~l (|2011D and 
IMorandi et al.l (|2011bf ) presented a triaxial model anal- 
ysis fitting simultaneously to the multiple-wavelength 
data of Abell 383, and Abell 1835, and derived the non- 
thermal pressure contribution, about 11% - 20%, out to 

-f?200- 

In this paper we continue our previous analysis of non- 
thermal pressure support in relaxed clusters assuming 
more general, non-spherical (ellipsoidal) geometry for 
the matter distribution. Since we are interested in the 
physics of the mass discrepancy, instead of comparing 
the total masses of clusters determined from X-ray and 
lensing observations, we test (as in M10) the assump- 
tion of strict hydrostatic equilibrium by comparing the 
two competing physical processes, self-gravity, acting in- 
ward, and the gradient of the pressure, acting outward. 
In particular, we study strict hydrostatic equilibrium in 
cluster MS 2137. 

MS 2137 (z = 0.313) has been studied exten- 
sively by using X- r ay and optica l (lensing) data 



(iHammer et all 119971: |G avazz 1 12001 IComerford et al] 
2006 1 : iShu et al.l 120081: iSand et al.l 120081: iMerten et al.l 
2009t IDonnarumma et al.l 120091: IComerford et al.ll2010D . 



Comerford et al. I (120061) modeled the central mass distri- 



bution of MS 2137 assuming an asymmetric NFW pro- 
file with the observed gravitational arcs. They found 
that the reconstructed lens system and the estimated 
mass profile are in a good agreement with simulations. 
IDonnarumma et al.l ()2009l ) (hereafter D9) estimated the 
mass of MS 2137 in the central and the outer regions by 
analyzing strong lensing data observed by Hubble Space 
Telescope and CHANDRA X-ray data, respectively. As- 
suming spherical and elliptical cluster model s with the 
major axis aligned in the LoS, iGava zzi 2005 (hereafter 
G5) derived the total mass profile in MS 2137 between 
0.02 and 1.0 Mpc by fitting the conventional and gen- 
eral NFW models to the strong and weak lensing data. 
G5 found a better agreement between masses determined 
from X-ray and lensing observations using prolate mod- 
els. G5 carried out a detailed combined weak and strong 
gravitational lensing analysis of MS 2137 (the same data 
were used by D9), therefore we use their results as our 
reference for lensing. 

MS 2137 is a relaxed cluster with a round morphol- 
ogy (close to circular) in projection (in the plane of the 
sky). In order to simplify the modeling, in this case, 
with a good approximation, we can derive the thermal 
and equilibrium pressure profiles for this cluster assum- 
ing spherical or the elliptical (prolate) models with the 
non-degenerated axis aligned in the LoS. 

The outline of our paper is as follows: in section 2, 
we discuss our X-ray data analysis of CHANDRA data 
for MS 2137, and the derivation of the thermal pressure, 
Pth , based on our analysis. We discuss our method to 
derive the equilibrium pressure, P eq , in section 3, based 
on published results from strong and weak gravitational 
lensing. We present our test of the assumption of hydro- 
static equilibrium in MS 2137, we discuss the possible 
sources of non-thermal pressure support and systematic 
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biases in section 4 and draw our conclusion in section 
5. Throughout in this paper, we assume a concordance 
cosmological model with il m = 0.3, ft\ = 0.7, and Hq = 
70 km s _1 Mpc -1 . Based on this model, the angular size 
1' = 0.275 Mpc for MS 2137. Unless otherwise stated, 
all errors, error bars and dashed lines arc represented la 
confident intervals. 

2. X-RAY DATA ANALYSIS 
2.1. X-ray Data Reduction and Cleaning 

We use publicly available observations of MS 2137 from 
the CHANDRA Data Archive (CDA). We found three 
observations with observation IDs (ObsID) 928, 4974 and 
5250 with exposure times of 44.17 ks, 58.14 ks and 41.09 
ks, respectively. We reduced these three data sets by us- 
ing the latest CHANDRA data analysis software package 
CIAO-4.4 and updated complement calibration CALDB- 
4.4. We followed the standard ACIS data preparation^ to 
reprocess all observations from level 1 event files (evtl) 
to level 2 reduced event files (evt2). First, in order to 
improve the quality of our data, we removed the streak 
by running the script destreak before we created the new 
bad pixel files. Then, we created the new bad pixel files 
correcting the known bad pixels and columns, identify- 
ing and marking the afterglow and hot pixels by running 
the scripts acis -find-afterglow and acis-build_badpix. Our 
next step was to run the script acis -process -events pro- 
ducing a new evtl with event grades 0,2,3,4 and 6 using 
the latest calibration data including the new bad pixel 
file created previously. Then we applied two stages of 
filtering to these new evtl files: first: for the imaging 
observation, we filtered for bad grades in the imaging 
observation, second: we used the primary Good Time 
Interval (GTI) to create the new evt2 file for each obser- 
vation data sets. In addition, we detected and removed 
the point sources by the script celldetect. Only a small 
number of the contaminated point sources was found re- 
sulting the removal of only a few pixels. 

Most of the contamination can be filtered out after the 
data reprocessing, but contamination from strong back- 
ground flares, which could increase the count rate several 
times, still has to be considered. Our final step was to re- 
move periods of high count rates due to background flares 
based on the light curves of each observation. However, 
the three observations have a different amount of contam- 
ination from strong flares. D9 studied MS 2137 by the 
joint-analysis of the arcs due to lensing observed by HST 
and CHANDRA X-ray data, but they only used ObsIDs 
928 and 5250, stating that there is a difficulty in iden- 
tifying the time intervals contaminated by background 
flares in the observation with ObsID 4974. Furthermore, 
the net exposure times in their work were reduced from 
44.17, 58.14 and 41.09 ks to 20.07, ~ 16 and 25.0 ks 
for ObsID 928, 4974 and 5250, respectively. As a result, 
they decided to discard ObsID 4974 in order to avoid 
the possibility of a systematic bias due to contamination 
from strong background flares. However, because the net 
exposure time contributed by ObsID 4974 is still con- 
siderable to the total net exposure time after cleaning, 
including ObsID 4974 can decrease the statistical uncer- 
tainty in our results. Thus, after a careful flare-cleaning 

2 http:/ /cxc. harvard.edu/ciao/guidcs/acis_data. html 



process, we decided to include ObsID 4974 in our study. 

The cleaning procedure for strong flare background for 
all three observations we adopted was the following. Af- 
ter the reprocessing the data creating evt2 files from the 
evtl files, the light curves of ObsID 5250 show no diffi- 
culty in filtering the background flares via running the 
CIAO-4.4 software script lc-dean. ObsID 928, on the 
other hand, is clearly contaminated by strong flares dur- 
ing the observation time about 59332 to 59352 ks, while 
the rest of exposure time shows a rather stable count 
rate. Therefore we discarded the time interval between 
59332 and 59352 ks for ObsID 928 before we run the 
script lc-dean. Parts of the S3 chip, 5' away from the 
emission peak, are selected as background fields in the 
script lc-dean for ObsIDs 928 and 5250, and the reduced 
net exposure times are consistent to the result of D9. In 
ObsID 4974 (with the net exposure time ~ 60 ks), we 
find strong contamination from background flares, which 
is also described in D9. In order to probe contamina- 
tion due to flares in ObsID 4974, we extracted the light 
curves from 5 annuli centered in the emission peak with 
radii uniformly distributed between 0.5' and 4.5' in the 
log space. We find that the light curves in these five re- 
gions show a similar trend: the count rate of the first 
~ 40 ks observation time is anomaly higher than adja- 
cent exposure time after that, and the count rate of the 
later exposure times, about 20 ks, is more steady imply- 
ing that this interval suffers from fewer flares. Based on 
our findings, we concluded that the background scatter- 
ing distributes approximately uniformly near the field of 
the cluster in the S3 chips but not regularly and periodi- 
cally spreading during the whole observation time. Other 
words, the background flares highly bias the cluster emis- 
sion in the first ~ 40 ks of the exposure, but much less 
after that. Hence, we decide to analyze only the last 20 
ks of the observation time for ObsID 4974. We further 
filtered the flares in this time interval using the region of 
the annulus, having the inner and outer radii of 2.89' and 
4.5', respectively, centered on the emission peak as our 
background criteria via running the script lc-dean. The 
GTIs for these three observations were determined after 
all corrections had been applied. The net exposure times 
after applying these GTIs are reduced from 44.17, 58.14 
and 41.09 ks to 14.07, 16.633 and 15.353 ks for ObsIDs 
928, 4974 and 5250, respectively. 

2.2. Background Subtraction 

The usual method to separate the X-ray emission due 
to galaxy clusters from the background is to subtract 
the background as determined from blank sky observa- 
tions from the data. However, CHANDRA observations 
have two telemetry modes, FAINT and Very FAINT 
(VFAINT) modes, which have different grading selec- 
tions. The standard blank sky background provided 
by the CHANDRA science team are reproduced from 
observations processed using the FAINT modes, while 
all three observations we are considering were using the 
VFAINT mode. Therefore, instead of making use of 
the blank sky backgrounds, we used local backgrounds 
with the assumption that the quiescent backgrounds 
did not vary significantly between the CCD chips. For 
each observations with ObsIDs 4974 and 5250, the back- 
ground was extracted individually from another back- 
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Fig. 1. — The dcprojcctcd temperature (upper panel) and the 
electron number density profile (lower panel) assuming spherical 
symmetry. The red color dots are the results of the simultane- 
ously fitting combining three observations and assuming spherical 
symmetry. The green curves are the best-fitted functions. 

illuminated chip (SI), which covers a part of the sky far 
from the cluster center, thus it can be considered free 
of cluster emission. Because of the lack of the SI chip 
in ObsID 928, we extracted the local background from 
S3 far away from the cluster center. As a consistency 
check, we compared this local background of ObsID 928 
to the background extracted from SI chip in ObsID 5250, 
which the data set suffers from less background flares. 
We found spectra of these two backgrounds are consis- 
tent with each other after normalization of the fluxes. 
Thus we apply this to background selection for ObsID 
928 and the local background of ObsID 928 is selected as 
the region of residual S3 chip having the distance of 4.5' 
far from emission peak. 



2.3. Spectral Fitting 

We extracted 9 circular annuli (8 complete and 1 frag- 
ment of the annulus having the spanning angle of 265° 
due to the lack of the field) centered in the X-ray emis- 
sion peak (2l' l 40 m 15.2 s ,-23°39'40") for all three data 
sets. We use the CIAO-4.4 software tool specextract to 
create the spectrum, the Response Matrix File (RMF) 
and the Ancillary Response File (ARF) for each region. 
Since MS 2137 is located at the center in one of the two 
back-illuminated (BI) chips, S3, in all three ObsIDs, we 
had no difficulty in extracting the spectra, RMFs and 
ARFs out to the radii of 220". The background spectra 
were extracted using the CIAO-4.4 software tool dmex- 
tract independently with the source spectra. 

We deproject the observed 2D spectral image to get 
the 3D gas properties. We adopt two types of models 
for the distribution for the ICG: spherical and ellipti- 
cal. Assuming circular symmetry in the plane of the sky 
and using a prolate model with the non-degenerated axis 
aligning in the LoS, we divide the gas distribution into 9 
shells according to the 9 annuli we used in the spectral 
extraction. Again, we assume the gas is homogeneous in 
each shell. 

We use the X-ray fitting package XS PEC- 12. 5 in our 
spectral analysis. We describe the gas emission in 
each shell by an optically thin plasma emission model, 



MEKAL, multiplied by the photo-electric absorption 
model, WABS. Our model for the projected spectra is 
generated by using the appropriate volume ratios associ- 
ating the normalization parameters in all MEKAL mod- 
els (see Appendix for more details of this projection). 
We fixed the galactic photoelectric absorptior0 and the 
rcdshift to the value 3.55 xl0 20 cm~ 2 and 0.313, respec- 
tively. We fitted for the three dimensional distribution of 
the temperature, abundance and normalization parame- 
ters using C-statistic in the energy band 0.5 to 7.0 keV. 

As a consistency test, first we fitted the spectra for 
each observation individually to see whether the results 
agree within errors. The temperature profiles for ObsID 
4974 and 5250 are consistent with the result of D9 except 
that our maximum values are slightly higher than theirs 
(we will discuss it in the next paragraph). On the other 
hand, the obtained deprojected temperatures and den- 
sities did not agree in the outer regions (above distance 
1.3' far from the emission peak) for ObsID 928, although 
the inner regions were consistent with the results derived 
from the other two ObsIDs and D9. We performed fit- 
tings to spectra in the outer regions (farther than 1.3') 
tentatively only in the 1.0 to 7.0 keV band (where the en- 
ergy bands of background spectra are the most consistent 
between the local and the blank sky backgrounds) using 
the blank sky background scaled to ObsID 928. How- 
ever, we found that there is no significant differences 
in the spectral fitting results between using the scaled 
blank sky background and the local background we se- 
lected (extracted from the region beyond 4.5') and the 
fitted energy ranges. The spectral fitting remained ex- 
tremely ill-constrained in the outer regions of ObsID 928 
possibly because the emission is dominated by the back- 
ground. In order to avoid bias in our final results, we 
decided to discard these regions (about 1.3' away from 
emission peak) of ObsID 928 when we performed simul- 
taneous fitting to all observations (i.e. we discarded the 
last three annuli in ObsID 928). 

We also compared the results of our X-ray data analy- 
sis to those of D9. There is a good agreement between the 
two density distributions (within errors). We compared 
our spectroscopic-like temperature profile projected from 
our best-fit 3D mode l using the weighting derived by 
iMazzotta et al.l (|2004h to the deprojected temperature 
profile derived directly from observations by D9 (sec red 
points in their Figure 3). Our temperature profile is also 
consistent with that of D9 except at about 200 kpc where 
our temperature is 1 keV higher than theirs. This dis- 
crepancy is most likely coming from two sources: i) We 
sample the outer regions (100-700 kpc) with more ra- 
dial bins than D9 (we are using 4, D9 are using 2 bins). 
At around 200 kpc there is a large temperature gradient 
which, if averaged over a large radial bin, results in a 
lower ring averaged temperature, ii) Statistical fluctua- 
tions. 

We found that there is no significant change in the dis- 
tribution of the temperature and density in our X-ray 
analysis for the assumed axis ratios for the dark mat- 
ter mass distribution, thus we only present our depro- 
jected temperatures and electron number densities for 
the spherical case (See Fig. [J). We use the virial radius, 

3 column density investigated from http://heasarc.gsfc.nasa 
gov/docs/tools.html 
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TABLE 1 

Best Fitted Temperature and Density Parameters 



To 
(kcV) 


9c 


Tc 

(Mpc) 


a c 


Tt 
(Mpc) 




11.6 
±1.37 


-0.65 
±0.0352 


0.107 
±0.009 


fixed to 
2.0 


0.478 
±0.079 


1.711 
±0.31 



no 

(10- 2 #/cm 3 ) 


n 

(lO" 2 Mpc) 


P 


9.24 ±0.43 


3.35 ±0.145 


0.61 ± 0.0074 



K v i r , as a distance unit in our figures to make it easier to 
compare our results to those of others (see the Appendix 
for the exact definition of our distance unit). 



2.4. Derivation of Thermal Pressure 

We use the following smooth parametrized functions 
for the dcprojcctcd temperature and density distribu- 
tions and perform fittings to the data points we obtained 
as described in the previous previous section: 

1 + Qr expj— (m„/r r ) ac } , . 
T 3D (m g (R,Z)) = T y % ■ Pt I 9 lj } } (1) 



n e (m g (R, Z)) = 



(l + m g /r T ) aT 
n 



(l + (m g /r b y)iP 



(2) 



where To is the temperature amplitude, g c is the scaling 
ratio, a c is the power index determining the falloff or 
increase of the temperature at the center, a? determines 
the falloff at the large radius, r c is the core radius and 
r<r is the temperature scale radius. The conventional 
beta model, Eqn. [31 is used to fit our density points. We 
define m g (R,Z) = H(to,€2,e3,x,y,z), as the elliptical 
coordinate, where 62 and £3 are the ellipticities of the 
gas, where R is the radial coordinate in the plane of sky, 
R = \/ y 2 + z 2 , and Z is the coordinate in the LoS, Z = z 
(note: £2 and 63 are functions of the radius). Because the 
results of the X-ray analysis are nearly insensitive to the 
ellipticities of the gas distribution we consider here, we 
show the best-fitted parameters of the temperature and 
density functions for the spherical case only (see Table. [TJ 
and we fixed a c = 2.0). The curves based on the best- 
fitted parameters assuming spherical symmetry for the 
temperature, electron number density profiles are shown 
in Fig. |U 

We use the best-fitted temperature and electron den- 
sity profiles to derive the distribution of the thermal pres- 
sure using the ideal gas law, P t h = nksTg, where P t h 
is thermal gas pressure, T g is the temperature of the gas 
assumed to be equal to the electron temperature and n is 
the total gas number density, the sum of the electron and 
ion densities, n = n e + rti (we assum^l n e /rii = 1.1737). 
The thermal pressure profiles for each mass model with 
different ellipticities arc shown on Fig. [5] 

3. EQUILIBRIUM PRESSURE 

In this section we describe our method to derive the 
equilibrium pressure, P eq , using published results from 
a combined analysis of weak and strong lensing observa- 
tions. The non-spherical geometry of the galaxy clusters 

4 see for example M10 



TABLE 2 

Best-fitted Mass Model Parameters 



q 


C200 


A/200 

(10+ 15 M Q ) 


vir 


vzr 

(10 +15 M Q ) 


(Mpc) 


1.0 


11.73 
±0.55 


0.772 

+0.047 
-0.042 


14.11 
±0.65 


0.847 
±0.0505 


2.05 


1.25 






12.75 
±0.52 


0.77 
±0.050 


1.85 


1.5 






11.60 
±0.50 


0.70 
±0.052 


1.68 


1.75 






10.744 
±0.477 


0.647 
±0.055 


1.556 


2.0 






9.92 
±0.406 


0.594 
±0.046 


1.446 



Note. — The mass model parameters. First column: the axis 
ratios q (axis aligned LoS over perpendicular to it) we choose 
to study, q = 1.0 implies the spherical symmetry. Second and 
Third column: the NFW model parameters A/200 and C200 in G5. 
Fourth, Fifth and Sixth column: converted spherical or elliptical 



NFW model parameters, 
virializcd states. 



and 



has been extensively studied by simulations or observa 
tions, and it has been suggested that the gas distribu- 
tion is most likely different to the tot al mass distribution 
assuming hydrostatic equilibrium dBuote fc Canizarel 
1998t " 



2004; 



mg n ydrost atic e quilibrium (Itsuote fe amzarcf 
■Ting fc Sutol \2(M ILee fc Sutol 1200$ IQguri et al 



Lau et al.1 1201 lh Therefore we derive the gravi 



tational potential using spherical and non-spherical ge- 
ometry assuming that the total mass density distribu- 
tion has fixed ellipticities, e2 and e 3 (which are not the 
same as the gas ellipticities: £2 = £2^2, e 3 , x, y, z) and 
£3 = £3(e2) &3, x, y, z)), and then obtain the gas density 
ellipticities and the gravitational force from the total po- 
tential. Finally the equilibrium pressure is derived as- 
suming the hydrostatic equilibrium equation. 

3.1. Spherical and Non-spherical Deprojection 

A gravitational lensing analysis provides the projected 
(2D) distribution of total mass surface density. In order 
to determine the 3D distribution of the cluster, we need 
to deproject the 2D mass surface density. 

For both cases, spherical and non-spherical symme- 
try, we use the more general elliptical NFW model 
parametrized by mAr(i?, Z) = LT(io, &2, 63, x, y, z): 



p N (m N (R,Z)) 



where po is the central density and r s is the scale radius, 
which is reduced to the spherical case for e 2 = e3 = 0. 

We use the published best-fitted NFW parameters 
( Af2oo ari d C200) from G5, and simulate the mass sur- 
face density maps based on their best-fitted parameters. 
Under our assumption of prolate models with the non- 
degenerated axis aligned in the LoS (i.e. e2 = 63), the 
surface mass density S is circularly symmetric in the 
plane of sky and can be expressed as the function of R 
as the following: 



E(iJ) = 



p N (m N (R,Z))dZ 



= 2 1 PN(m N ) 

lm=R 



qm N 



-.dm 



N 



(4) 
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-H-| 1 1 1 1 — I I I I | 1 1 1 PN' I I I 




M-| 1 1 1 1 — I I I I | 1 1 1 1 t" l I I 




1.00 



Fig. 2. — Pressure profiles as a function of projected radius (in 
the plane of the sky) for different axis ratios q of the dark matter 
mass distribution (from top to bottom: q = 1.0, 1.25, 1.5, 1.75 
and 2.0). We show the observed pressure distribution (red points 
with error bars), and the best fitted profiles of the thermal and 
equilibrium pressures (green and blue solid and dashed lines). 



q is the axis ratio (axis in the LoS over perpendicular 
to it) defined by q = 1/ y/l — e% = 1/ y/l — e| for our 
prolate models. For q = 1.0(ore2 = P3 = 0), Eqn. [4] is 
reduced to the spherical geometry. 

For our prolate models, we estimate the elliptical NFW 
model parameters - the elliptical virial mass M^ ir and 
the elliptical concentration parameter G e vir - by fitting 
the models to the mass surface density map. We found 
that there is no significant effect on M^ ir and C e uir 
from how we choose the simulated data points for fitting, 
thus we generated 15 simulated S data points distributed 
uniformly from 0.01 to 3.0 Mpc for our elliptical fitting. 
The errors for the best-fitted parameters, M^ ir and C^ ir 
, were generated from 1000 Monte-Carlo realizations. We 
show our best-fitted mass model parameters in Table. [2j 
Once M^ ir and G e vir are obtained, we define the elliptical 
virial radius R^ ir based on the average density enclosed 
within the prolate ellipsoid with projected radius R^ ir , 
as 

M^ T (R e Mr ) 3 g Pm (z)A(z), (5) 

where p m (z) is the mean matter density and A(z) is the 
virial overdensity at redshift z. 

3.2. Derivation of Equilibrium Pressure 

We derive the equilibrium pressure P eq using the equa- 
tion of the strict hydrostatic equilibrium: 



(6) 



where p g is the gas mass density and $ is the total gravi- 
tational potential. We derive the gravitational force from 
the total gravitational potential $. For a given NFW 
mass distribution pm and a set of the total mass el- 
lipticities e-i and e^, $ ca n be written as the following 
(jBinnev fc Tremaind[200l l: 



(7) 



Jo 



00 
t=0 

n 2 



dt 



tp{oo) - v>(n 2 ) 



y/T+Ty/t + 1 - + 1 



dm 2 N p N (m 2 N ) 



For our case of prolate models, Eqn. [6] can be param- 
eterized using R, and reduces to 



1 dP, 



(8) 



p g dR 

where p g is the gas mass density. Using Eqn. g(R) can 
be written as the following: 

4irG f qR m 2 N pN{rnN)drriN 

' mjM—O 



I 



-4irGpor s q x 

q R 



Ri + mlil-jz) 



xdx 



{l+xf (^y +x 2 {q 2 _ X) 



(9) 



where x = m^/r s and G is Newton's constant. For 
q = 1, Eqn. [9] reduces to the expression for spherical 
geometry. 
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Fig. 3. — Pressure ratio profiles Pth/Peq for MS 2f37 assuming 
axis ratios of q = 1.0, 1.25, 1.5, 1.75 and 2.0 (red, green, blue, pur- 
ple and black solid and dashed lines). The horizontal grey dashed 
line indicates hydrostatic equilibrium. 

Using the best-fitted NFW model parameters and the 
gas density profiles derived from X-ray analysis (assum- 
ing the mean molecular weight fj, = 0.61879, as in M10), 
the equilibrium pressure, P eq , profiles can be derived 
numerically by integrating Eqn. [8j Nevertheless, Eqn. |8j 
is a first-order differential equation thus one boundary 
condition is needed for a solution. Unfortunately, the 
pressure is not known neither at the cluster center, nor 
at the virial radius (where it is not decreasing to zero 
due to the accretion shock) . 

Analyzing AMR and SPH simulations, Molnar et al. 
have found that the accretion shocks in clusters are lo- 
cated at 0.9-1.3 R V i r , where the pressur es in different 
clust ers are a few times 10~ 13 dyne/cm 2 ()Molnar et al.l 
l2009t 2010). Currently, we cannot derive the pressure 
at K v i r due to the very weak X-ray signal from large 
radii, therefore we have no constraints on any boundary 
conditions from observations. Fortunately, the pressure 
at R vir , P eq (R = Rvir), contributes less than 0.1% to 
the central pressure, therefore the choice of its value has 
no significant effect on the pressure in the inner part 
of the cluster. However, deriving the pressure distri- 
bution based on different boundary conditions (different 
P eq {R = Rvir)), we noticed that the slope of the pressure 
profile changes only close to R V i r , and it is about about 
0.8 near R V i r = 0.7. Therefore assuming the same slope 
for radii close to R V i r provides a smooth pressure dis- 
tribution, which we would expect out to R V i r (see M10 
for the example). The equilibrium pressure profiles we 
obtained based on the assumption above are shown in 

Fig. m 

4. DISCUSSION 

4.1. Testing Hydrostatic Equilibrium 

We test the assumption of strict hydrostatic equilib- 
rium in MS 2137 assuming spherical and prolate models 
by analyzing the pressure ratio profiles R pr = Pth/Peq- 
We use five different axis ratios for the dark matter dis- 
tribution (the semi-axis along LoS over the perpendic- 
ular to it): q = 1.0 (spherical case), 1.25, 1.5, 1.75, 
and 2.0. R pr = 1.0 implies strict hydrostatic equilib- 
rium, while R pr < 1 implies that non-thermal pressure, 



Fig. 4. — Non-thermal pressure ratio profiles, Rnpr = 1 — 
Pth/Peq = Pnth/ Peq, for MS 2137 with different axis ratios (the 
assumed axis ratios and color codes are the same as in Fig. [3] We 
also show Rnpr profiles of A1689 and two simulated relaxed clus- 
ters studied in MlO under the assumption of spherical geometry 
(orange and cyan solid and dashed lines). 

P n th = Peq — Pth, is required for the cluster to be in 
equilibrium. R pr > 1 would imply that the thermal 
pressure is larger than the pressure required for equi- 
librium, therefore the cluster is expanding, which is not 
acceptable for relaxed clusters (there would be no phys- 
ical reason for the cluster to expand). We show the Rpr 
profiles as a function of projected radius R in Fig. |3l 
All ratio profiles arc plotted only within the radial range 
where both X-ray and lensing observations were available 
(no extrapolations were used). 

Assuming a spherical model for MS 2137, the thermal 
pressure contribution is about 50% ± 5% within 0.03 
R v ir , and it increases to a maximum of about 60% at 
0.08 Rvir , and decreases again to about 25%-40% at 
0.2-0.5 R V ir ■ In MS 2137, overall, the thermal pressure 
contribution is almost half of the equilibrium pressure at 
all radii, which is consistent with the result of G5. As a 
result, the mass determined from X-ray analysis by as- 
suming spherical symmetry and hydrostatic equilibrium 
is systematically less than the mass based on the lensing 
analysis by 40% at all radii. As it was shown by G5, 
a possible reason for this discrepancy is non-spherical 
geometry. We find that, assuming prolate models with 
different axis ratios, q (principal axis aligned with the 
LoS), the non-thermal pressure contribution necessary 
for equilibrium is reduced, and thus the two mass esti- 
mates would be closer. From Fig. [31 q — 1-25 can result 
less non-thermal contribution at all radii and be clos- 
est to the hydrostatic equilibrium at about 0.1 R V i r , 
implying that MS 2137 could possibly be far from the 
spherical symmetry Assuming q = 1.75, we can re- 
solve the discrepancy within 0.03 R V i r for MS 2137 but 
non-physically overestimating the P t h at 0.07-0.2 R V i r 
. Setting q = 2.0 for MS 2137 overestimates the Pth 
significantly at almost the whole radii. However, for the 
large radii beyond 0.5 R V i r , MS 2137 is deviated signifi- 
cantly from the hydrostatic equilibrium among q ranging 
from 1.0 to 2.0. 

We show the non-thermal pressure ratio profiles, 
Rnpr = I- Rpr (or P nth = P eq - P«,), for MS 2137 
in Figure ([4]). For comparison, in this figure, we also 
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show non-thermal pressure ratio profiles of Abell 1689 
and simulated clusters using results from M10 (where 
we assumed spherical symmetry). The R npr values for 
MS 2137, assuming spherical symmetry, are systemati- 
cally higher than the those derived from simulations by 
2 er, providing further evidence that spherical symmetry 
for MS 2137 is not a good assumption. Assuming q = 
1.25-1.5, R np r of MS 2137 we find a good match with 
the simulations within 0.1 R OT > , but they deviate at the 
large radii (beyond 0.15 R V i r ), suggesting that the non- 
thermal pressure components could possibly contribute 
significantly at the large radii for MS 2137. Although the 
assumption of hydrostatic equilibrium in MS 2137 seems 
to be better if we adopt q = 1.25-1.5 (statistic value 
is 11740.0 with 10632 degrees of freedom, reduced \ 2 = 
1.10 for q = 1.25), we still find a significant non-thermal 
contribution of about 10%-40% in the core region. Our 
results suggest that, even if we assume elliptical mod- 
els, strict hydrostatic equilibrium is not valid in the core 
region of galaxy clusters, and therefore the contribution 
from non-thermal pressure has to be taken into account 
when accurate cluster mass determination from X-ray 
observations is preformed. 

Interestingly, the R n pr profile for A1689 seems to 
be qualitatively different from the profile derived for 
MS 2137 and from those derived from simulations. While 
the R np r profiles for MS 2137 and those from simula- 
tions monotonically decrease to their minimum at about 
0.1 R-uir , and increase at large radii, R npr for A1689 is 
roughly constant in the center region (within about 0.06 
K V i r ), decreases to its minimum at about 0.25 R„i r , 
and increases monotonically for large radii. The min- 
imum non-thermal pressure contribution in Abell 1689 
located at around 0.25 R V i r is about twice as far as those 
in MS 2137 and simulated relaxed clusters. Our analysis 
of pressure profile ratios suggest that Abell 1689 is not 
a typical relaxed cluster. Possible explanations for the 
peculiar pressure ratio distribution in Abell 1689 can be 
due to substructure/ the lack of dynamical equilibrium 
in the core of Abell 1689, and/or feedback from heating, 
which could push the most relaxed regions outward in 
the cluster. Clearly, more detailed studies are needed to 
understand the structure of Abell 1689. 

4.2. Sources of Systematic Bias 

Since we are comparing the equilibrium pressure, 
which depends directly on the lensing measurement, 
to the thermal pressure which depends directly on the 
X-ray data, our results arc sensitive to the calibra- 
tion of the two measurements. In practice, we find 
that most of our systematic errors arc coming form 
the lensing mass measurement. Several sources of sys- 
tematic bias in lensing measurements have been iden- 
tified, such as the inherent shapes of the background 
galaxies, signal dilution due to cluster galaxies, the 
low signal-to-noise ratio, substructure, contamination by 
mass concentrations along the LoS, aspher i city, the as- 
sume d mass models (IRozo fc Schmidtll2010HOguri et al.l 
20101: iQguri fc Takadal 120111 : iSchmidt fe Rozd 12011 



Assuming spherical clusters, the dominant system- 
atic bias is originated form triaxiality, which may re- 
sult an err or in the mass measurement more than 
40 - 50% (ICorless fc Kind 120081: [Hamana et al.l I20T2I: 
iSereno fc Zitrinl 120121 IBahe et al.l 120121 ). In order to 
avoid this main source of uncertainty, we determine pres- 
sure profiles using ellipsoidal models with different val- 
ues for the ellipticity, q , and to have a conservative 
estimate on the non-thermal pressure support, we ac- 
cept the value of q , which provides the least amount of 
non-thermal pressure support without violating our as- 
sumption of equilibrium. This method is not sensitive to 
systematic bias in the total mass determination due to 
substructure or mass concentrations in the LoS. 

On the other hand, the X-ray analysis is not affected 
significantly by the assumed geometry of the cluster. 
However, clumping of the gas in the outskirts would en- 
hance the X-ray emission resulting in ov erestimate in the 
temp e rature and the dens i ty of the gas (iMathiesen et al.l 
199| INagai fc LavJ l20ll lEckert et al.l I2012D . Shocks 
would result an excess in the X-ray emission and an 
overestimate of the thermal pressure, resulting a lower 
non-thermal contribution. Of course, our assumption of 
equilibrium would be wrong for merging clusters. We 
choose to use MS 2137 specifically to avoid most of these 
systematic biases. MS 2137 shows no sign of a recent 
merger event, no significant substructure, and it looks 
nearly circularly symmetric in projection. We avoid sys- 
tematic bias due to extrapolation of our parametrized 
models to radii we do not have data by studying the 
cluster in the radial region observations cover. Therefore 
we do not expect significant bias in our results. 

4.3. Sources of Non-Thermal Pressure Support 

Several physical mechanisms have been proposed for 
non-thermal pressure support in the last few decades. 
The most likely sources are magnetic, cosmic ray pres- 
sure and pressure due to turbulent gas motion. Mag- 
netic fields in the clusters could grow from the hi- 
erarchical formations or could be injected by galactic 
winds, the n possibly further undergo galactic dynamo 
proce sses dGoncalves fc Friacal ll999TlDolag fc Schindlerl 
l2lMlGwoni fc Ferettill2004t IChurazov et al.ll2008T) ~ 

Magnetic fields in the core of several clusters of 
galaxies have been estimated using Faraday rota- 
tion measurements, and found to be in the range 
of 10-30 (2-8)u G fo r (non-)cool core c lusters (i.e., 
Tavlor fc PerTevl Il99l iFeretti et all H991 iTavlor et al.1 



2002f ). It has been suggested that t he magnetic field 



Gruen et al.ll2011t iBecker fc Kravtsovll20lH IBahe et al.1 



2012f h The accumulated systematic errors due to the un- 
certainty of the background galaxies and the mass con - 
centration outside -R200 is about 10% (|Bahe et al.ll2012l ). 
resulting the same order bias to the pressure ratios. 



strength is proportional to t he de n sity dDolag et al.lll99£ 
iColafrancesco fc Giordanol l2007t lAndo fc Nagail l2008Tl . 
Using simulations, iDolag fc Schindlerl (|200dT T have found 
that the averaged magnetic pressure is about 5% of the 
thermal pressure in relaxed clusters, implying that other 
sources of pressure are needed to explain the discrepancy 
between the X-ray and lensing mass estimates. 

Shocks generated by structure formation, active galac- 
tic nuclei (AGNs) and supernova feedback may accel- 
erate particles to very high energies. It has been sug- 
gested that cosmic rays generated by these processes 
could also provide a source for non-thermal pressure sup - 
port (lEnsslin et al1H997t [Saraiin1[200l [&oofcOl][200a . 
Constraints on the cosmic ray pressure contributions can 
be placed by studying the gamma rays emission gener- 
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ated by neutral pions produ ced by collisions betw een cos- 
mic ray and ICG particles (|Ando fc N agai 2008;). Based 
on numerical simulations and observations, the contri- 
bution from cosmic rays pressure was estimated to be 
less than 10-30% of the thermal pressure in the very 
central region (with i n about few tens k p c) of clusters 
(lEnsslin et all 119971: iMiniati et all 120011: iReimer et all 



20031: IPfrommer fe Ensslinl 120041: IPfrommer et all 120071 
200 a IPfrommerll2008l: ISiiackT et al. [2008D. 

The most widely accepted candidate for the non- 
thermal pressure support is turbulent gas motion. Cos- 
mological simulations show that even relaxed clus- 
ters contain 100 kpc scale subsonic flo ws and turbu- 
lence generated b y stru ct ure formation (lEvrar dlH990l: 



Norman fe Bryan 
20071 : IFang etal 



INagai. Vikhlinin fc Kravtsov 
ILau et al.l 120091 : iVazza et al 



20111 ). At the outskirts of clusters subsonic flows are 



most likely generated by previous major mergers, and ad- 
ditional turbulence can be generated by supersonic mo- 
tion of cluster galaxies. Accretion via filaments delivers 
gas and galaxies to the core region of clusters providing 
power for turbulence in the core. Non-gravitational pro- 
cesses such as jets and bubbles from AG Ns may also con- 
tribu te to turbulence in the core region (jChurazov et all 
120021) . Based on simulated clusters, M10 showed that 
even without additional non-gravitational physics the 
non-thermal pressure support in the core of relaxed clus- 
ters can be about 30%, reaches a minimum in the core re- 
gion (at about 0.1 R V i r ), and increases to about 35% at 
the virial radius. M10 also showed, that assuming spheri- 
cal models, there is a significant contribution, about 35%, 
from non-thermal pressure at the central region of A1689. 
In this paper we have shown that, even if we use elliptical 
cluster models, the strict hydrostatic equilibrium is not 
valid in MS 2137, we find a significant, 30%-40% contri- 
bution to non-thermal pressure in the core region. Our 
result indicate that the significant non-thermal pressure 
support in the core region of clusters is not a consequence 
of the assumed spherical shape of the clusters, but a gen- 
eral feature. 

A detailed study of the sources of different non-thermal 
pressure contribution in MS 2137 would require more de- 
tailed observations in different wavelengths, and there- 
fore it is out of the scope of our paper. In prin- 
ciple, the different non-thermal pressure contributions 
could be estimated using parametrized functions, but 
the results would str ongly depend on the parametrization 
(jLagana et al.l 120101 ). The quality of MS 2137 observa- 
tions available to us would not allow a meaningful separa- 
tion of the non-thermal pressure contributions, therefore 
we have presented our results in terms of thermal and 
non-thermal pressure support. 

5. CONCLUSION 

We have derived thermal and non-thermal pressure ra- 
tio profiles for cluster MS 2137 from X-ray and gravita- 
tional lensing observations assuming spherical and ellipti- 
cal models. Since MS 2137 show nearly circular projected 
morphology, with a good approximation, in our analysis, 
we could assume elliptical (prolate) models with the non- 
degenerate axis aligned in the LoS. 

Analyzing the pressure ratio profiles, we have found, 
in accord with previous studies based on different argu- 
ments, that a prolate model with q = 1.25-1.5 for MS 



2137 provides an adequate description for the distribu- 
tion of the ICG. We have found that the contribution 
from non-thermal pressure varies with radius. Based on 
our analysis of the pressure ratios ( R pr and R npr ), 
spherical symmetry does not seem to be a good assump- 
tion for MS 2137. Assuming q = 1.25-1.5, R npr reaches 
its minimum value of 0%-20% at about 0.08 R vir , and 
the non-thermal component contributes about 10%-40% 
and 20%-60% at small (within 0.03 R V i r ) and large 
(0.2-0.5 R V i r ) radii, respectively. We conclude that 
MS 2137 is not in hydrostatic equilibrium in these two 
regions even assuming non-spherical symmetry. Our re- 
sults from comparing the shapes of pressure ratio profiles, 
Rnpr , for MS 2137, A1689, and those derived from re- 
laxed clusters drawn from cosmological simulations, sug- 
gest that A1689 is not a typical relaxed cluster. 

Our results based on X-ray and gravitational lensing 
observations, and numerical simulations of relaxed mas- 
sive galaxy clusters provide further evidence that there 
is a significant contribution from non-thermal pressure 
in the core region of even relaxed clusters, i.e., the as- 
sumption of hydrostatic equilibrium is not valid in this 
region, independently of the assumed shape of the cluster 
(spherical or elliptical). 

INC thanks K. Umetsu for discussions regarding the 
lensing analysis, and H.-H. Chan for advice on deprojec- 
tion techniques. We thank the referee for comments and 
suggestions which helped us to improve our paper. 
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APPENDIX 

We project our gas models onto the plane of sky assuming that the gas is in the equilibrium with the gravitational 
potential generated by our assumed total mass distribution. In our mass model, the dark matter density distribution 
is described by a series of mass shells with constant axis ratios, q (principal axis along LoS over the minor axis, 
which is perpendicular to LoS). The equipotential surfaces generated from each mass shells coincides with a series of 
the confocal ellipsoidal surfaces, and the axis ratio for these equipotential surfaces is getting closer to 1 (the shapes 
are becoming rounder as the radius increases). We derive the total potential by summing up the potential generated 
from each mass shell. As a consistency test, we checked how much the total potential changes once we add the gas 
to the system, and found that the change is insignificant, less than 0.1%. We concluded that the geometry of the 
equipotential surfaces is dominated by the dark matter halo, so we can ignore the effect of the gas. We assume the 
gas distribution is identical to shape of the total potential by the X-ray Shape Theorem (jBuote fc Canizareslll998f ). 

Before projection, we define our coordinate system which is used throughout this paper. We define the confocal 
ellipsoidal coordinates, n(t, e 2 , e 3 , x, y, z), by the following: 

II 2 (i, e 2 , e 3 , x, y, z) = JL- + ^f— | + TT ^— | (1) 

e| = l - (a 2 /a 1 ) 2 
e 3 = l - (a 3 /ai) 2 , 

where a\,a 2 and a 3 are the three semi-axes and a 3 < a 2 < ai (i.e. < e 2 < e 3 < 1). x, y and y are the coordinates of 
three semi-axes, t is the parameter representing the family members of the same series of the confocal surfaces. Given 
a fixed e 2 and e 3 , the shape of II surface becomes rounder as t increases. For e 2 — e 3 — 0, II is reduced to the spherical 
distribution. For t = constant, II describes a conventional ellipsoidal distribution. 

The emissivity in clusters is dominated by the thermal bremsstrahlung radiation and is equal to the product of the 
square of the gas number density, n 2 , and the cooling function, A. For convenience, we label the inner most to the 
outer most annuli, where we extract the spectra from, by 1 to N, consecutively. Assuming the gas is homogeneous 
within each shell, the expected luminosity we observed for the fcth annulus in the plane of sky, L k , can be expressed 
as: 

L k = Y.tunelKVk,* , (2) 

where i represents the quantity of the ith shell, and Vi j is the volume of the jth shell occupying the cylinder having the 
non-degenerated axis aligned LoS and the base of the ith annulus. By introducing the volume ratio Rk,i = Vk^/Vi^, 
the Eqn. [5] can be reduced to the following form: 

L k = '£g k (n e 2 i A i V i , i )R k , i (3) 

Eqn. [3]can be viewed as the basic ideal of the spectral fitting. Vij as well as Rij = 1 , 2, . . . , N) can be computed 
directly from the gas distribution. We derive our gas models in the 3D space and project them onto the plane of sky for 
a given set of volume ratios Rk,i, which connect the normalization parameters in the MEKAL models of the different 
shells. We fit the models to the spectra simultaneously to get the temperature, the abundance and the normalization 
parameter of each shell. 
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